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Abstract 

We consider a Kondo-like impurity interacting with fermions on a honeycomb lattice at half- 
filhng, as in the case of graphene. We derive from the lattice model an effective one-dimensional 
continuum theory which has, in general, four flavors with angular momentum mixing in the presence 
of internode scattering processes and six couplings in the spin-isotropic case. Under particular 
conditions, however, it can be reduced to a single-coupling multichannel pseudogap Kondo model. 
We finally calculate, in the presence of an energy dependent Fermi velocity, induced by Coulomb 
interaction, the critical coupling in the large-A^ expansion, the magnetic susceptibility and the 
specific heat. 

PACS numbers: 71.10.Fd; 72.15.Qm; 75.30.Hx; 71.10.Ay 
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I. INTRODUCTION 

Since the experimental realization of a single monolayer of graphite named graphene, 
a two dimensional crystal made of carbon atoms hexagonally packed, a lot of efforts have 
been made to study many properties of electrons sitting on a honeycomb lattice {2 1 . Also the 
problem of magnetic impurities in such a system has become a topic of recent investigations 
in the last few years , although a detailed derivation of the effective model is 

still lacking. 

The main motivation of the present work is, therefore, that of deriving, from the lattice 
Hamiltonian, the corresponding continuum model for the Kondo-like impurity, writing the 
effective couplings from the lattice parameters. From angular mode expansion we get an 
effective one-dimensional Kondo model which has, in general, four flavors and is peculiar 
to graphene-like sublattice systems. Strikingly, we find that there is an angular momentum 
mixing only in the presence of internode scattering processes, being the valleys and the 
momenta locked in pairs, in each sublattice sector. The complete model has six couplings 
in the spin-isotropic case, however, thanks to the lattice symmetry, for some particular 
positions of the impurity, the number of couplings can be reduced to one, obtaining a 
multichannel pseudoga p K ondo model sharing, now, many similarities with other gapless 

'ermionic systems js], yJlO], as for example, some semiconductors 8|, d-wave superconductors 

11 1 and flux phases 12|. 



A second issue which is worthwhile being addressed is related to interactions. In real 
systems logarithmic corrections in the density of states may appears, as a result of many- 
body effects. In order to include, at some extent, correlation effects we allow the Fermi 
velocity to be energy dependent. Indeed for a system of electrons in the half- filled honeycomb 
lattice, like graphene, an effect of Coulomb interaction is that of renormalizing the Fermi 
velocity |13i] which grows in the infrared limit. This behavior induces in the density of 
states subleading logarithmic corrections. We plan therefore to analyze the effect of these 
corrections onto the Kondo effect in order to see how finite coupling constant transition, 
obtained within the large-A^ expansion technique jl4] and renormalization group approach 



161], can be affected by deviations from power law. We find that the critical Kondo 



coupling becomes non-universal and is enhanced in the ultraviolet by a quantity directly 
related to the Coulomb screening. Moreover, we find that the impurity contribution to the 
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magnetic susceptibility and the specific heat vanish faster by log^ than in the free 
approaching zero magnetic field or zero temperature. 



II. THE MODEL 

In this section we will derive the continuum one- dimensional effective model from the 
microscopic lattice Hamiltonian. 

A. Lattice Hamiltonian 

Let us consider a honeycomb lattice which can be divided into two sublattices, A and B. 
The tight-binding vectors can be chosen as follows 

S,^^{1,V3), (1) 
S,^^{l,-V3), (2) 
53 = a(-l,0), (3) 

where a is the smallest distance between two sites. These vectors link sites belonging to two 
different triangular sublattices. Each sublattice is defined by linear combinations of other 
two vectors, |(3, \/3) and |(3, —VS). From these values one can derive the reciprocal-lattice 
vectors in momentum space and draw the Brillouin zone which has an hexagonal shape, i.e. 
with six corners. We choose two inequivalent corners (the others are obtained by a shift of 
a reciprocal-lattice vector) at the positions 

K' = -K. (5) 

These points are actually the Fermi surface reduced to two dots approaching the zero chem- 
ical potential, i.e. at half-filling. 

We will consider the following Hamiltonian defined on this honeycomb lattice 

H = Ho + Hk. (6) 
The first contribution is given by the tight-binding Hamiltonian 

Ho^-t^c\^{v)cBa{Y + 5)+h.c., (7) 

rScr 



where t is the nearest neighbour hopping parameter, c\^{r) (cyio-(r)) is the creation (anni- 
hilation) operator for electrons with spin cr localized on the site r, a vector belonging to 
the sublattice A, while c^^^{r + 6) (cBcr(r + 5)) the creation (annihilation) operator for elec- 
trons on the site r + 6, belonging to the sublattice B. The second contribution to H is the 
Kondo-like impurity term 

Hk = Y1 {^v± (5+c{(v)q(v) + S.c\{v)c^{v)) + K^S, {c\{y)c^{y) - c\{y)c^{v)) } , (8) 

V 

where Av± and Av^ are the short-range Kondo couplings, S (with S± = Sx ± ^5"^) is the spin 
of the impurity sitting at the reference position (0, 0), a is the spin operator of the electrons 
located at v from the impurity, v can belong to A or i? and we sum over all these vectors. 

B. Derivation of ID effective model 

We now rewrite the fields c in the following way 

c^.(r)^e^^->],,(r) + e-^^->L(r), (9) 
CB.(r + 5) ^ e*^-(^+^VL(r + S) + e-^^-('-+^VL(r + S), (10) 
c.(v) ^ e^^->L(v) + e-^->L(v), (11) 

with L = A if V G A, or L = i? if V G 5. The upper indices, 1 and 2, label the Fermi points 
Eqs. (jl]), ([5]). At these particular points we get the following equalities 

^e±'^-^ = 0, (12) 

<5 

E^^""' = -T(1'T0. (13) 

<5 

Expanding the slow fields 'ip\^^{r + 6) around r, introducing the multispinor 

the identities do, tq, 70 and the Pauli matrices cxj, and 7^, i = 1,2,3, acting respectively 
on the spin space, t|, valley space, 1, 2, and sublattice space. A, B, we get, in the continuum 



limit, 



Ho = -ivp j (ir^'^(r)o-o {r-ilidy - ro729^) V'(r), (15) 
Hk = ^^(0) Qjx (^+a_ + S_a+) + J.S,a}j ^(0), (16) 



where vp = ^ is the Fermi velocity and 

4 = ^ { {Joiro + Jtri + ^21^2) (70 + 73) + ( J(f.ro + jf^n + J^r,) (70 - 73) } (17) 

a Kondo couphng matrix with the following components, containing the lattice details, 

4. = E^v., (18) 
Ji^, = 5^cos(2K-v)Av., (19) 
J2^, = ^sin(2K-v)Av<,, (20) 

vS-L 

where a =_L, z the spin index and L = A,B the sublattice index. Eq. (fT5l) is a Dirac- 
Weyl Hamiltonian, constant in spin-space, which, after defining ip = ip'^'y^, can be written 
as Vp J driplr) {rQ'-fidx + T-^'~f2dy) ipir) to make Lorentz invariance manifest. The spectrum 
is made of a couple of Dirac cones departing from the two Fermi points, and the density 
of states vanishes linearly approaching the zero energy, p(e) = z/|e|, where u oc Vp'^. This 
property plays a fundamental role on the scaling behavior of the Kondo impurity, as we are 
going to see. Let us rewrite the full effective Hamiltonian in momentum space, 

^^-p ^"^(p) p (To (r37i sin 9p - ro72 cos 9p) z^(p) 
+ / / ^^q) (S^-- + S-a^) + IS.a.^ ^(p), (21) 

where we have parametrized the momenta as follows 

= p cos 9p, (22) 
Py = p sin 9p. (23) 

For the benefits of forthcoming discussions we first notice that the orbital angular momentum 
operator 

C = ~t{xdy - yd,) (24) 

does not commute with the Hamiltonian Hq in Eq. f|T5l) . On the other hand, in order to 
define proper total angular momenta we introduce the operator 

J = /: + ^r373, (25) 
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which does commute with Hq, 



[J,Ho]=0, 



(26) 



and also with the ro-components of Hk- In particular, given some amplitudes '4'Yam{p)i 
eigenstate of J with eigenvalue j = (m + can be written as 



Performing the following unitary transformation 



(27) 



U„ 



2v^ 



■To [(1 + te-''^^'){^o - ^72) + (1 - ze-^^-^«)(7i - 73)] 



to the fields 



^(p) = Up(f){p), 



the Hamiltonian Eq. ( l2Tl) becomes 



H = vf 
+ 



dp f (iq 



PJPcro^o73</>lPJ 



(28) 



(29) 



(30) 



'>t(q) ( ^K^{9^, Op) + + K,{9^, ^p)5,(T, ) 0(p), 



(27r)2 ; (27r)2 

namely, i^o becomes diagonal, the cost to pay is that the Kondo couplings depends on the 
angular part of the momenta, 

+ (jo^.et(^-^^)-ro + jf,e-t(^^+^^)-n + JZe-'^^'^^'^^^r,) (70 + 7i)}.(31) 

Notice that the angular dependence of K does not prevent the model to be renormalizable. 

16| . being q the intermediate momen- 



As one can see by poor man's scaling procedure 
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tum in the edge bands, dropping for the moment the spin indices, the contributions which 



renormalize, for instance, Jq in the particle channel are 



2 i^Pl ^V2 



)f3[ 



(32) 
(33) 
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with i = 1,2. In the same way we can check that the corrections to Jf, with i 

Analogous corrections can be verified in the hole channel. In all these corrections 6q always 
cancels out, recovering the right momentum dependence for the slow modes. 
From Eqs. fll8H20l) we actually get access to the renormalization of linear combinations of 
the original lattice parameters Av. 

In order to reduce the problem to one dimension we proceed expanding the fields <f){p) 
in angular momentum eigenmodes as follows 

CO 

0(p) = J2 e'^'^^-^^^'^Mp), (36) 

m=co 

with m G Z. Indeed, due to the gauge in Eq. (!28l) . all the spinor components have the same 
angular phase. Actually from Eq. ( l29l) . one verify that e*''™'"*"5)^p0„(p) is the eigenvector of 
J', Eq. (1271) . with eigenvalue j = m + |, transformed by Up^, and with amplitudes 

AamiP) = ^ ii^BamiP) T , (37) 

where the subscript ± replaces the sublattice index and refers to the sign of the energy, 
vfPIs, appearing in Eq. (130!) . The original field at position r = (r, ip), can be written as 

^(^ ^) = / ^ y ^rn^^rn^ ( ' (P^e^^T-i + Uprhu] 0^(p) \ ^^^^ 

J ^V27c^t^^ y lUprhd- lm+i{pr)e'''lu](pUp) j ' 

where ^rn{z) are the Bessel functions of the first kind, 7^ = 70 + 71 — ^72 — 73 and 7^ = 
7o ~ 7i ~ "^72 + 73- At r = the only terms which survive are those with m = 0,-1, 
corresponding to j = ±|, in terms of eigenvalues of JT. After integrating Eq. fl30|) over 
the angles, indeed, we get in the Hk only contributions with m = 0,-1, in the following 



= 1, 2, are 

(34) 
(35) 
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combinations 
H = 



dp 



VfP 



1=1.2 m=—OQ 



dp 



dq 



<Pl\q){S ■ a){^o - iMip) + • a)(7o - 7i)0'ib) 

+Ji [4HQ)iS ■ ^)(7o - 7i)0' i(p) + <P'M<1){S ■ a)(7o - 7i)0o(p) 
-^4 [<l>l\q){S ■ a){^o - 7i)0-i(p) - <P'Mq){S ■ a)ho - iMip) 

+Jo [<l>'Mq)is ■ a)(7o + hHUp) + 4Hq)is ■ ^){io + hHUp) 

+Jf WMq)iS ■ a)(7o + 7i)0o(p) + <PoiQ)iS ■ <?)(7o + 7i)0-i(p) 



— iJ. 



(39) 



where now (plnip) are spinors only in spin and energy spaces. Here we are considering the 
spin-isotropic case, with Jj" = J^j_ = J^^, to simphfy the notation. In the spin- anisotropic 
case one simply has to replace ■ a) with Jli_{Sx(Jx + Syay) + J^^^SzCTz)- In the free 

part of the effective model, Hq, we keep only the contributions from the particles with 
m = 0, — 1, the only ones which can scatter with the impurity. We now unfold the momenta 
from [0, +oo) to (— oo, +oo) by redefining the fields in the following way 



Csaip) = [sign(p)]' 



\p\(Plign{p)am{\p\), 



(40) 



where, in order to label the fermions, we choose the index i in valley space and the index 
s = sign (m + I) , the sign of the total angular momenta, eigenvalues of J', which are good 
quantum numbers as soon as there is not internode scattering, i.e. = = J2 = J2 = 0. 
Introducing for simplicity 



j2Kv 



A. 



(41) 



we finally end up with the following one-dimensional effective Hamiltonian 

s- {4 {ejiqWM+f-^i^inup)) + j^ejiqweip) + j^e(gX(p) 
+ Jo (tHQ)^e-{p)+&{<i)^e4p)) + j^tHqWM + JU'KmHp) 



(42) 



where, in the first term, the indices s = ±, i = 1,2 and a =1, i, are summed, and the 
dispersion relation is E{p) = vpp. The full model, Eq. (j42l) . has six Kondo couplings, in the 
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spin-isotropic case, which are independent for a generic position of the magnetic impurity 
on the lattice. Moreover Eq. fH2|) exhibits an angular momentum mixing in the presence of 
internode scattering amplitudes and J^, namely, when also the nodes are mixed. We 
are not going to analyze the complete model in full generality but we shall consider only 
particular cases physically relevant. 

C. Some particular examples 

a. Impurity on a site. If we consider an impurity on top of a site of the honeycomb 
lattice, belonging to the sublattice A, for instance, and consider only the nearest neighbour 
coupling between the impurity and the electrons located on this site, we have Av 7^ if 
V = (0, 0) and assume Av = for v 7^ (0, 0). In this case we get 

Jq = Ji = A(o,0), (43) 
J2 = Jo = J? = Ji = 0. (44) 

Introducing the symmetric combination for the fields 

c = a + e!, (45) 

the effective Hamiltonian Eq. fH2l) becomes simply 

H = J2^E{p)Y,Cl{p)Up) + f //3|^iv^^-CngR(p), (46) 

which is a single channel Kondo model. 

b. Impurity by substitution. If we now consider an impurity sitting on a site of the 

honeycomb lattice, let us say, belonging to the sublattice A, and consider only nearest 

neighbour couplings between the impurity and the electrons, we have Av = if v G A while 

Av = A^j = A52 = A^a, if V = 6i, i = 1,2, 3, and Av = for t> > a. Noticing that 

3 3 

cos(2K ■6i) = J2 sin(2K ■ 6,) = 0, (47) 

i=l i=l 

we get a remarkable reduction of the number of couplings given by Eqs. fll8H20l) 

Jo^ = J^ = J^ = Jf = Jf = 0. (48) 



9 



Recalling the fields as follows 



(49) 



the effective Hamiltonian Eq. fj42l) reduces to 

Nf 



i=l 



where the Nf = 2 flavors (the valleys and the momenta are locked in pairs) are decoupled 
and we realize a two-channel Kondo model. The reduced model Eq. (l50l) is the same as that 
found for flux phases 121]. 

c. Impurity at the center of the cell. Finally, let us consider an impurity at the center 
of the honeycomb cell. In this case, using Eqs. fll8H20l) and Eq. P7|) . we have 



jA _ jB 



(51) 
(52) 



Enumerating the fields as follows, for instance. 



Ci — ^+ , C2 — ^- , Cs — , C4 — C 



2 

+ ' 



(53) 



we get the same Hamiltonian as in Eq. ( l50l) with, now, Nf = 4 flavors, realizing, therefore, 
a four-channel Kondo model, as in the case of (i-wave superconductors [ll|. 



III. LARGE- EXPANSION AND THE ROLE OF COULOMB INTERACTION 



In this section we solve the model Eq. fl50|l in the large- iV approximation, where N is 
the rank of the symmetry group of the impurity, which actually is equal to 2 for spin one- 



half. Following the standard procedure ll|, llj]; within a path integral formalism, we write 
S = flaa/sfiB, introducing additional fermionic fields /, with the constraint Q = flfa, the 
charge occupancy at the impurity site. In the Lagrangian, therefore, a Lagrange multiplier 
eo is included to enforce such constraint, which is actually the impurity Fermi level. To 
decouple the quartic fermionic term one introduces the Hubbard- Stratonovich fields 
where i = 1, ■■■Nf, being Nf the number of flavors. For impurity by substitution Nf = 2, 
as seen before. After integrating over the fermionic fields, ( and /, we end up with the 
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following effective free energy, 

^=7/ defie)5{e) + J dr | |$,(r)p - Q eo J , (54) 



where /(e) is the Fermi function and 



, n\e\A/2 \ . , 

d[e) = arctan -^5 . , , , ,. ., (55) 

the phase shift, with A = - |$j(e)p/7r, and A a positive ultraviolet cut-off which dictates 
the limit of validity of the continuum Dirac-like model for the free Hamiltonian. For graphene 
the typical value is A ~ 2eV. 

So far we have considered a model of free fermions hopping on a lattice and scattering 
eventually with a magnetic impurity, but in order to get more realistic predictions we should 
consider, at some extent, interaction effects. In order to do that, we let the Fermi velocity 
be energy dependent, i.e. vp = vpie). 

This is not unrealistic since it has been shown [13] that, due to Coulomb screening in an 
electronic system defined on the half-filled honeycomb lattice, as in the case of a monolayer of 



graphene 
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18| , the effective Fermi velocity is renormalized in such a way that vp flows to 
higher values in the infrared, and consequently the density of states around the Fermi energy 
decreases. The low energy behavior for the renormalized velocity is vp ~ ln(e~^) and so the 
density of states should behave naively as p ~ evp'^ ~ e/ ln(e^^)^. The aim of the following 
section is then to study the role of such corrections onto the Kondo effect, neglecting, 
however, possible renormalization of the Kondo coupling due to Coulomb interaction. The 
idea is to consider an uncharged magnetic impurity embedded in a cloud of charges dressed 



by Coulomb interaction. The realistic expression for the Fermi velocity is the following 

^^(6)=^(l + r^ln(A/|6|)), (56) 

where 77 is related to the flne structure constant, for Thomas- Fermi screening it is 77 = 
e^/Aehv, being e the dielectric constant, and v is an energy independent velocity. 

A. Saddle point equations 

From Eq. fl54p . the extremal values of 69 and A, evaluated at zero energy in the static 
approximation, satisfy the saddle point equations |^ = and |£ = 0, which can be written 

11 



as follows 111, ll2| 



TT 



Q = 1 rrf6/(6)^(6), (57) 

^ def{e)^{e), (58) 

" J-D 



where D < A is the bandwidth. The Eq. fl57j) dictates the relation between the singlet 
amplitude A ~ ^.-(l Xlo-Ctr/o-P) ^"^^ the impurity level eo, at fixed occupation charge Q, 
and reads 

^ TT ^(7r|e|A)2 + 4(.;^(e)2(e-eo) + eAln(A/|e|))2- ^ ' 

For T = and for a generic value of Q, we get the following behavior for the impurity 

A(l + >7ln(A/g)) 



level, Co ~ A e'' ^^^2»7Qu^(i+i7in(A/D))-Ay _ jj^ ^j^e non-interacting limit, formally, when 77 — 0, 
it reduces to eo ~ D e"^"^*^/^, in agreement with Ref. 11]. Strikingly, the limit of A 
is finite and equal to Ae^/^, i.e. the two limits do not commute. This means that, in the 
presence of Coulomb interaction, the occupation charge for an impurity level within the 
bandwidth is finite only if the singlet is formed and Q ~ A/2'r]vvF{D). The energy scale 
eo in Eq. (!59|) does not play the role of an infrared cut-off for A — 0, and as a result, in 
that limit, Q goes to zero for any value of eo. This result is different from that found in the 
free case [11|] where the Fermi velocity is constant, VF{e) = v. In the latter case eo vanishes, 
approaching zero singlet amplitude, for any value of the occupation charge. 
The second equation, Eq. (1551) . dropping the indices for simplicity, reads 

1 ^de^.,, 2t;^(e)2|e|(eo-e) 



J 7r^^'^(7r|e|A)2 + 4(t;^(e)2(e-eo) + eAln(A/|e|))2- 

Setting eo = and A = at T = 0, we get the following critical coupling 

1= r ^^_ = ^L_ [ --e^/''r[0,l/r/ + ln(A/Z})]l, (61) 

where r[a, x] = t"'~^e~^dt is the Incomplete Gamma function. Sending 17 — > we recover 
the standard result -j- = js]. 

At this point it is worthwhile making a digression. Contrary to the free case, where the 
limit lim^i^o ^fj- is trivially finite and equal to in the interacting case, using Eq. (jSI]), 
this limit is zero. On the other hand, if we replace v with renormalized velocity vf{D), the 
limit 

lim ^ = ± (62) 
12 



is finite and equal to the standard case. This is consistent with the fact that the dimension- 
less parameter relevant in the Kondo effect is not the bare coupling J but the product pJ 
and that, in the presence of a renormalized Fermi velocity, Eq. fl56p . the density of states 
is modified as p(e) ~ e/fp'(e)^. In order to validate this result and to get more insights one 
can address the problem from a renormalization group prospective, as we did in Appendix. 
From Eq. (l6T!l we find that the critical coupling Jcp{D) is not universal, being an increasing 
function of the ratio D/A, and is larger than the corresponding mean field result in the 
non-interacting case for any positive D < A. 

To go beyond the tree level, one should consider quantum fluctuations, i.e. higher orders in 
large- expansion, which might spoil the critical point obtained in the mean field level, as in 
the case of strictly power-law pseudogap Kondo systems pl)|], if the particle-hole symmetry 
is preserved. In order to break particle-hole symmetry, however, one can include straightfor- 
wardly a gate voltage in the model In any case the role of fluctuations, in the presence 
of a logarithmic deviation from power-law in the density of states is still an open issue which 
we are not going to address here. 



B. Magnetic susceptibility and specific heat 

d. Magnetic susceptibility. The magnetic field can be easily included in our final model 
introducing a Zeeman term Ha^. This term modifies the phase shift in the free energy, 
Eq. (15^ . as 5(e) — > |(5(e + H) + 6{e — H)). We can, therefore, calculate the magnetization 

M{T,H) = -^, (63) 

and the magnetic susceptibility 

x(T,i/) = -^. (64) 

For T ^ and if ^ eo, we have the following magnetization 

Af AAA AAA 

M(0, H) ^ - IH-H) - Sm - IMHr ^AMA/H)] . ■ 

(65) 

The final result for the magnetization is valid only ii H Ae~^^^'"'^^'^ . In the same limit the 
asymptotic behavior of the magnetic susceptibility is, then, given by 

^^^^""^^ elwHmr ^^^^ 
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For Ae-^/(''^)' < < eo, instead, one gets x(0, H)^ e[vfiH)^ H^/H), and for vf{H) 



111 



one recover the result for the non-interacting case 
For if — > and T <^ cq, Ae^^/^^'^-'^, we have, instead, the following magnetic susceptibility 



ae^e 2e2 (1 + e-)3 t;^(T|x|)2 e^v'^r]^ ln(A/T)2' 

(67) 

which crosses over to x{T,0) = ^|^^^T ln(A/T), for Ae^^/^'''')' < T < eo- In the 
presence of Coulomb interaction, therefore, the magnetization and the susceptibility vanish 
logarithmically faster, as approaching zero magnetic field or zero temperature, than in the 
free pseudogap case. The overscreening effects pointed out in Ref. is, then, enhanced 
by a factor on the order [In ( ^^^^^ h) ) ] ^ ' presence of a renormalized Fermi velocity 

induced by the Coulomb screening. 

e. Specific heat. Let us calculate now the impurity contribution to the specific heat, 
defined as follows 



CiT,H) = -T^. (68) 



For H and T ^ cq, we have 



N ^^df d6 _ NAT^ ^ f x'^e''{e'' - 1) 2xe^ ^^ \x\x 



C(T,0) = — / dee'^—^—^ dx 



TnJ_^ dede 2el J \ [l + e-f [l + e-f ) VF{T\x\y 

9C(3)A^A T2 



,(69) 



elv^ ln(A/T)2' 

where C(3) ~ 1.2 is the Riemann zeta function at 3. Also in this case we assume T <^ 

9C(3)jVA- 



Ae ^/(''^)^. For T larger than that energy scale, instead, C(T, 0) ~ ^^^^^^^T2ln(A/T), and 



for ?7 ^ one recover the non-interacting result 



For T ^ and H <^ we have, instead, the following behavior 



7r2 



yTx(0,i7), (70) 

with x(0, H) calculated before. Eq. ( |70l) corresponds also to the asyptotic behavior of the 
impurity entropy, —dF/dT, for T <^ H. The specific heat, like the magnetic susceptibility, 
vanishes logarithmically faster than in the case with constant Fermi velocity, approaching 
zero magnetic field and zero temperature. However, the Wilson ratios C/{Tx), both for 
T <^ H and for H <^T, are exactly the same as those found in Ref. jllt |. 
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IV. SUMMARY AND CONCLUSIONS 



We have derived the low-energy continuum hmit of a Kondo-hke impurity model defined 
on a honeycomb lattice at half-filling. 

By angular momentum eigenmode expansion we have obtained an effective one- 
dimensional model with two colors and four flavors, two for each sublattice sector, Eq. (l42l) . 
The impurity effective Hamiltonian involves two angular momenta which are linked with 
the two nodes. We have found, therefore, that the internode scattering contributions cor- 
respond also to the angular momentum mixing terms. Quite in general, we have to deal 
with six couplings in the spin-isotropic case, which are linear combinations of the original 
lattice parameters, Eqs. fll8ti20p . However, due to the underlying lattice symmetry, in tight- 
binding approximation, the number of Kondo couplings can be reduced to one, for particular 
impurity conflgurations. 

We have flnally calculate, both within large- expansion technique and renormalization 
group approach the mean fleld critical Kondo coupling which is increased by the presence 
of a renormalized Fermi velocity driven by Coulomb interaction. From the calculation of 
some thermodynamic quantities, we flnd, however, that even though the Kondo phase is 
suppressed, at least in the mean fleld level, once the singlet is formed, Kondo screening 
effects are enforced by the Coulomb charge screening. 
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APPENDIX A: COMPARISON WITH RENORMALIZATION GROUP 

Redeflning the Kondo coupling as J' = with u a constant factor, we are going to 
consider the following generalized Kondo equation, 

dJ' dlnp I ,2 



de de 



J' + J'\ (Al) 



where £ = \n{D/t) with e > 0, a positive deflned energy parameter and p the density of 
states, not yet deflned. We shall be seeing that, quite in general, — lim^^oo^^ is the 
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infrared point at criticality which may or may not be a fixed point. Solving Eq. (lAip and 
requiring that J{D) = Jo we obtain 

J'{e) = -^'"P^'l , . ■ (A2) 

For positive couphng, i.e. antiferromagnetic Kondo model, we can define a critical point as 

and rewriting Eq. (lA2p in terms of we have 

J'(e) = J'oJcPi^) (A4) 

From these result one can immediately see that, if < then J'(e) for e while if 
J'o > J'^ then J'(e) oo at some low energy scale, called Kondo scale Tk- In particular Tk 
(let us put the Boltzmann constant fc^ = 1) is defined, for J'^ > J'^, by 







o c 



At J'o = J'^ the infrared limit is given by 



lim J\e] 



= = lini = - hm (A6) 

^ () X 



10 X 

where we have used the de I'Hopital rule in the second equality. 

Notice that we get really an unstable infrared fixed point from Eq. ( ]A4[) only for a special 
functional form of the density, let us call it p, which is the solution of 

Pie) = J'J^ (AT) 
Jo X 

which is p(e) = z/e"^^, i.e. a power law. 

As an example, let us consider logarithmic deviation from power law of the following form 



ln(^ 



(A8) 



e 

such that for r = a = we should get the standard Kondo result and for r 7^ and a = 
we should recover the known results for gapless systems. Assuming the form given by Eq. 
(IA8I) for the density of states, and using Eq. flASP we get the following result for the critical 
value 

J^ = r"+i±^ L^^^^ J (A9) 

A-r [l + a,rln(A)]' ^ ^ 
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which is not universal, dependent (increasingly for r > and a < 0) on the ratio D/A. 
From Eq. flA4l) we have the following expression for the coupling 

r"+i(^^-^^)Pp) + ^«A^r[l + a,rln(A)]- ' 

For smaller than J'^ it flows to zero, while for J'^ greater than J'^ it diverges approaching 
the following Kondo scale, calculated from Eq. (jA5|, 

( T - Wl"^ ( r \r-,, ji ji -\^ -"A 



K 



(All) 



From Eq. (lAlOp and using the fact that lim^^o a^j — ^ A^r~" or, equivalently, using 

Eq. flA6p for p given by Eq. flASp . we obtain that the infrared limit right at the critical 
point (i.e. at = J^) is lime^o = ^i which is the same as without logarithmic 

corrections. We see that this result does not match with given by Eq. (]A9p . except for 
a = 0, since, in this case, F [l, r In (^)] = jr? and we recover the known result Jc = r which 
is then an unstable fixed point. 



Motivated by a physically relevant example 
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18j, we now use a more realistic density 



of states induced by the renormalized Fermi velocity Eq. (1560 . which then reads 

P(e) = 2 ( 1 + tttttI • (A12) 

(l + r7ln(A/e)fV l + vH^/^)J 

Notice that in the infrared limit Eq. ( ]A12p reduces to Eq. ( lASp with r = 1, a = —2 and 
u —>■ v/rf. Using Eq. flA12p in Eq. flASP we get the following non- universal critical Kondo 
coupling 

" r]VF{D)p{D) -vue^/'^AT[0,l/r] + \og{A/D)y ^ ' 

which is almost equal to Eq. ( !6TI) , after rescaling. It becomes exactly the same if we approx- 
imate Eq. (1A12P with p(e) = J^'^^^^^^a , i-e. neglecting the second term. 

For 7] — ^> 0, in the non-interacting limit, we recover = 1. Also in the infrared limit 
D we get J'^= 1. Moreover, for D going to zero Eq. flAlSP approaches Eq. flA9p with 
r = 1 and a = —2, which does not depends on 77. On the other hand, for Z) ^ A we have 
•^c = nW^-l^^r[o,i/r,] ' ^nd for small r] we get 

~ 1 + 27]. (A14) 
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Namely, for D < A, in order to enter the Kondo regime, one needs a starting coupling 
which exceeds the critical coupling, obtained in the non-interacting case, by approximately 
an amount directly related to the Thomas- Fermi screening length. 
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